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Abstract 

A few issues on gluon confinement are addressed with the help of a renormalizable gauge model 
obtained by introducing a replica of the Faddeev-Popov action and a soft breaking of the BRST 
symmetry. Confinement turns out to be encoded in the spectral properties of the corresponding 
correlation functions. While the propagators of the elementary fields have no interpretation in terms 
of physical excitations, examples of local composite operators whose two-point correlation functions 
possess a spectral representation with positive spectral density can be introduced. These composite 
operators turn out to be left invariant by the BRST transformations, a feature which has strong 
consequences on their renormalizability properties. Moreover, they display a direct interpretation in 
terms of i-particles [TJ , which are the unphysical modes corresponding to a confining propagator of the 
Gribov type. A possible way to take into account the effects of the Gribov copies is also outlined. 
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1 Introduction 



Color confinement is a fascinating topic, representing a great challenge from theoretical, numerical and 
experimental points of view. Although impressive results have been achieved in the last decades, it is safe 
to state that many aspects remain still to be unraveled, forcing us to provide new ideas and to improve 
our current computational capabilities. 

The aim of the present work is twofold. Our first goal is that of investigating confinement within the 
framework of a local and renormalizable Euclidean quantum field theory. Here, confinement is meant to 
be encoded in the analyticity properties of the correlation functions, as expressed by the Kallen-Lehmann 
spectral representation, which provides a powerful tool in order to attach a physical meaning to the cor- 
relation functions [2j[3]. More precisely, in a confining theory, we expect that the two-point correlation 
functions of the elementary fields, here generically denoted by {(/>}, cannot be cast in the form of a spectral 
representation with positive spectral density. A meaningful particle interpretation of the propagators of 
the theory cannot thus be achieved. As a consequence, the excitations corresponding to the elementary 
fields {(/>} do not belong to the physical spectrum of the theory. As a familiar example of this, let us 
quote the Gribov propagator [4], which we shall employ extensively throughout: 

(WW-*)) = . 

where k 2 > is the Euclidean momentum and 7 is the so called Gribov mass parameter. This kind 
of propagator arises in the Gribov-Zwanziger theory [H [5j [6j [8] , which implements in a renormalizable 
way the restriction of the domain of integration in the Euclidean Faddeev-Popov quantization formula 
to the interior of the region in field space bounded by the Gribov horizon^. Expression (pQ) displays 
complex poles, so that it cannot be interpreted as the propagator of a physical particle, being suitable for 
a confining phase. It can be seen as describing the propagation of two unphysical modes with imaginary 
masses ±i7 2 , which have been called i-particles in [T], namely 

1 

+ 



fc 4 + 7 4 2 \k 2 -ij 2 k 2 + ij 2 

Furthermore, a suitable set of local composite operators, {0(<p)}, should be introduced in such a way that 
the two-point correlation functions, (0(k)0(—k)) , exhibit the Kallen-Lehmann spectral representation 

(0(k)0(-k)) = j™ dr , (3) 

where p(r) is the positive spectral density and tq > stands for the threshold. As it follows from complex 
analysis, expression allows us to define an analytic function in the complex cut £; 2 -plane [2J[3], where 
the cut extends along the negative real axis, from — tq to —00. An important outcome of expression 
([3D is represented by the welcome feature that we can move from Euclidean to Minkowski space, i.e. 
k 2 — > —k 2 M , where the cut will be now located along the positive real axis, starting at the threshold +to 
and extending till +00. Moreover, positivity of the spectral density p(r) enables us to give an interpreta- 
tion of eq.([3]) in terms of physical states with positive norirH. The operators {0((p)} would thus provide 



1 See ref.[9] for a general account on the Gribov-Zwanziger theory. 
2 The ideal situation would be that in which, in addfl 
a real pole in the correlation function (0(k)0{ — k)) , i.e 



2 The ideal situation would be that in which, in addition to cuts along the negative real axis, we would be able to isolate 
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information on the spectrum of the theory. 

One should notice that the correlation function (0(k)0(—k)) is supposed to be evaluated by employing 
a confining propagator of the Gribov type, eq.(pQ), which exhibits complex poles. It is thus not obvious 
that expression ([3]) will display the desired analyticity properties. 

To have an idea of the difficulty related to the use of a confining propagator of the Gribov type, let 
us quote the case of the two-point correlation function (F 2 (x)F 2 (y)}, F 2 (x) = F^ v (x)F®(x), which has 
been investigated at one- loop order in [6] [7]. In spite of the fact that F 2 (x) is a BRST invariant and renor- 
malizable operator within the Gribov-Zwanziger theory it gives rise to a correlation function which, 
besides exhibiting a cut along the negative real axis, does have unphysical cuts along the imaginary axis 
[6] . A promising way out in order to avoid the presence of the unphysical cuts has been outlined in [1] , 
where the introduction of the i-particles has been proven to be helpful in order to construct examples of 
composite operators whose correlation functions exhibit the Kallen-Lehmann spectral representation at 
one-loop order. Though, it turns out that, unlike F 2 (x), these composite operators are not left invariant 
by the softly broken BRST symmetry which ensures the renormalizability of the Gribov-Zwanziger action 



Further work seems to be needed in order to achieve a deeper understanding of the relationship between 
the softly broken BRST symmetry and the construction of a sensible set of renormalizable composite 
operators exhibiting good analyticity properties. This is the second topic we aim to investigate in the 
present paper. We shall elaborate on a model which, in analogy with the Gribov-Zwanziger theory, dis- 
plays a softly broken BRST symmetry. Moreover, we shall be able to introduce examples of composite 
operators which are invariant under the softly broken BRST symmetry, while having a direct interpreta- 
tion in terms of i-particles. This feature enables us to immediately establish the Kallen-Lehmann spectral 
representation for the corresponding correlation functions at one-loop order. In addition, the invariance 
property of these operators translates into softly broken Slavnov- Taylor identities which play a central 
role in order to establish how these operators renormalize, i.e. how they mix under radiative corrections. 
The knowledge of the mixing matrix is needed in order to correctly evaluate the correlation functions 
at higher orders as well as to find out quantum operators which are invariant under the renormalization 
group equations, see [12] [11] and references therein. 

We point out that the issue of the existence of a BRST symmetry compatible with the confining character 
of the theory is one of the highlight of the current investigation on the infrared behavior of Yang-Mills 
theories and, in particular, on the possible existence of a non-perturbative version of the BRST symmetry. 
We remind here that this symmetry plays a fundamental role in the proof of the renormalizability and of 
the unitarity of non-Abelian gauge theories in the perturbative regime. At present, it is unknown if and 
how such a powerful symmetry can be reconciled with the existence of the Gribov copies. What emerges 
from the studies of the Gribov-Zwanziger theory is that the BRST symmetry enjoyed by the Faddeev- 

where M is a mass parameter and 1Z(k 2 ) an analytic function. Such an expression would provide a direct link with the 
physical spectrum of the theory. The mass M would refer to the mass of a physical excitation created by the action of the 
composite operator 0((j>) on the vacuum of the theory. This is what one expects for the spectrum of the glueballs in QCD, 
which should manifest as poles in the correlation functions of suitable gauge invariant composite operators [10]. Nevertheless, 
from an analytic point of view, finding a pole is a highly non-trivial matter. To the best of our knowledge, this would require 
a re-summation of an infinite set of diagrams [21 [3J, a not easy task indeed. A more modest plan, which will be assumed 
here, would be that of not facing the too hard calculation of the pole, limiting ourselves to the determination of the cuts 
along the negative real axis, something which can be already regarded as a non-trivial achievement, given the complexity of 
the Gribov propagator Q. Though, we might also argue that the knowledge of the threshold to could provide us a rough 
estimate of the mass M, as follows by observing that to should be related to the energy for the creation of multi-particle 
states [2] [3j. 
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Popov action is softly broken by terms which are proportional to the Gribov parameter 7 |13LI14| H3]. As 
such, the breaking of the BRST symmetry seems to be unavoidably linked to the presence of the Gribov 
horizon. Despite the various efforts done recently [Ml H3 [THJ, [H], the issue of the BRST versus the 
Gribov horizon is certainly far from being understood. Therefore, it seems to us that the construction 
of a confining model exhibiting a softly broken BRST symmetry and allowing for the introduction of 
examples of composite operators which are BRST invariant and whose correlation functions display the 
Kallen-Lehmann spectral representation might be of some help for the current debate. 

The paper is organized as follows. In Sect[2]we present the model and its softly broken BRST symmetry. 
We motivate its introduction by showing that it has a natural interpretation in terms of i-particles [I], 
i.e. in terms of the unphysical modes corresponding to the confining Gribov propagator, eq.([T]). In Sect[3] 
we provide an example of a local composite operator whose two-point correlation function exhibits the 
Kallen-Lehmann representation at one-loop order. This composite operator turns out to be left invariant 
by the softly broken BRST symmetry, a feature which is relevant for the study of its renormalizabil- 
ity, which we shall briefly outline. SectJH deals with the issue of the Gribov copies. A possible set up, 
relying on the use of Gribov's no-pole condition [JJ, is analysed in order to account for the Gribov phe- 
nomenon. In Sectj6]we present a refined version of the model which, in analogy to the so-called refined 
Gribov- Zwanziger action introduced in [20, [13], gives rise to a modified gluon propagator attaining a finite 
non-vanishing value at zero momentum, a behavior also reported by the analysis of the Schwinger-Dyson 
equations [211 E21 [231 02] as well as by models employing an effective gluon mass [2H [25] . We mention 
here that recent lattice numerical simulations [261 123 (2H1 [291 EOl EU [32] point towards a gluon propagator 
exhibiting a finite non- vanishing value at zero momentum. Sect [7] collects our conclusion. A few technical 
details about the Ward identities fulfilled by the model are given in the Appendix [Al 



2 The model and its softly broken BRST symmetry 

The construction of the gauge model which we shall present in the following relies on the mastering of 
the confining Gribov type propagator, eq.([T]), which we have been able to acquire so far [I]. 

2.1 A confining gauge model 

Let us start by considering the Faddeev-Popov action in the Landau gauge 

S FP = J d 4 x (±F° U (A)F^(A) + b a d,A; + c a d^ b (A)c b ^j , (5) 

where b a is the Lagrange multiplier enforcing the Landau gauge condition, d^A® = 0, (c a ,c a ) stand for 
the Faddeev-Popov ghosts, F^ V (A) is the field strength 

F; U {A) = d,Al - d v Al + gf^A^Al , (6) 

and D^(A) denotes the covariant derivative 

DfiAy^d^ + gr'A^ . (7) 

Our confining model is obtained by performing the following steps: 
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• we first consider a replica of the Faddeev-Popov action, eq.Q, by introducing a set of mirror fields 
(Z7° ,b a , u) a , uj a ) as well as a mirror Faddeev-Popov action: 

Smfp = J d 4 x (^U^(U)U^(U) + b a d^ + u a d^Df(U)u b ^j , (8) 

where U^ U (U) is the field strength corresponding to the field U® 

KAU) = fyttf - d u U; + gf abc U b U c u , (9) 

while b a and (ij a ,uj a ) stand for the mirror Lagrange multiplier and mirror Faddeev-Popov ghosts, 
respectively, and 

D ab {U)J = d^ a + gr b U c ^ h . (10) 

• the gauge field A^ is softly coupled to its mirror field U£ through the following mixed term 

S e = iV20 2 J d 4 x A^U® , (11) 

where 6 is a mass parameter which, for the time being, is considered as a free parameter introduced 
by hand. As we shall see, this mass parameter will play a role akin to that of the Gribov parameter 
7 of the Gribov-Zwanziger action [H El [8] . Later on, in SectJH we shall attempt to provide a 
possible dynamical mechanism for 9. 

Our confining model is thus specified by the following action 
S = Sfp + Smfp + Se 

= J d'x (\f^(a)f^(a) + \u; u {u)u; u {u) + ^^x) 

+ j d 4 x ^Al + ^d^DfiA^ + ^d^ + ^d^DfiU^y (12) 
A first property of the action (|12p follows by looking at the propagators of the fields (A®, U%), i.e 

(A a ,(k)AU-k)) = ^^2F " kj w) ' (13) 

K(k)U b u (-k)) = S^j^; (V - ^) , (14) 

(A«Wt(-k)) = (<W - ^) • (15) 

As one sees from expressions (fT5|) . JHJ), (US]), all propagators are of the confining Gribov type. As such, 
they correspond to unphysical excitations. Said otherwise, we cannot attach a particle interpretation to 
the propagators of the elementary fields (A 1 ^, U®). 

The second feature displayed by the model is that it contains a unique coupling constant g. Both 
fields A® and £/" interact with the same coupling. We shall provide further comments on this aspect 
when the relationship with the z-particles will be discussed. Nevertheless, we point out that the feature 
of having a unique coupling constant is protected by a powerful discrete invariance which we shall refer to 
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as the mirror symmetry. It turns out in fact that the action S, eq. (ll2p . is left invariant by the following 
discrete transformations 



Ad 
\ ~ 


v TT a 






u;- 




b a - 


+ b a 


i a - 


-> b a 






c a - 


> u a 






oo a - 


-> c a 






(f - 




uj a - 





(16) 



The mirror symmetry (116p means that the fields (A^, b a , c a , c a ) can be replaced by (U^,b a ,O a ,u; a ), and 
vice- versa. Notice that the mirror symmetry would be lost if a second coupling, g' , would have been 
introduced. Suppose in fact that the field strength of the field U£ would contain a different coupling, i.e. 

U« V (U) = d.U? - d v u; + g'f abc U b U c u . (17) 

Evidently, the mirror symmetry would not hold anymore, as U£ V (U) cannot be transformed into the 
mirror field strength F®(A). 

A third feature displayed by the action (112p is the existence of a softly broken BRST symmetry. It 
is easy to verify that the nilpotent BRST transformations 

sAl = -D ab (A)c b , 
sU; = -Df(U)u\ 

sc a = ^-f abc c b c c , 
2 

a = 9jabc b c 

2 

s(f = b a , 
sb a = 0, 
su a = b a , 

sb a = , (18) 
leave the action S invariant up to soft terms proportional to the parameter 9 2 , i.e. 

sS = 9 2 A break , (19) 

where A break is given by 

A break = -iV2 J d*x (u«D ab {A)c b + AlD ab {U)j) . (20) 
Being of dimension two in the fields, A break is a soft breaking. 
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2.2 Recovering Yang-Mills theory in the deep perturbative ultraviolet region 



An important aspect of the Gribov-Zwanziger action is that the perturbative ultraviolet behavior of 
Yang-Mills theory is recovered when the Gribov horizon is removed. More precisely, when the Gribov 
parameter is set to zero, 7 = 0, the Gribov-Zwanziger action can be proven to be equivalent to the 
Faddeev-Popov action. This property of the Gribov-Zwanziger action is encoded in the BRST trans- 
formations of the auxiliary fields (<^f, <^f, ui^ , 0^ b ) needed to implement the restriction to the Gribov 
region [9j. These fields form a BRST quartet, thus yielding a trivial BRST cohomology. As the BRST 
symmetry becomes an exact symmetry when 7 = 0, it follows that the BRST quartet decouples from the 
theory, so that the Faddeev-Popov action is recovered at the end, see [S]. 



This property is also displayed by our action 5, eq. (|12j) . although the recovering mechanism is deeply 
different from that of the Gribov-Zwanziger theory. One notices that, when 8 = 0, the BRST transforma- 
tions of the mirror fields (U^,b a ,Q a ,uj a ), eqs. (j!8j) . remain non-trivial, i.e. they give rise to nonvanishing 
BRST cohomology classes [33]. However, it turns out that, when 9 = 0, the mirror fields (U^,b a ,u) a ,uj a ) 
do not interact with (A't, b a , c a , c a ). In fact, by construction, when the mass parameter 8 is set to zero, the 
two set of fields (A®, b a , cf, c a ) and (Z7*, b a , u) a , u a ) belong to two disjoint theories which do not interact, 
namely 

S = S FP (A,b,c,c) + S M Fp{U,b,uj,uj) . (21) 

As a consequence, denoting by {0\{A, b, c, c), On(A, b, c, c)) a set of operators depending only on the 
original fields (A" b a , c°, c a ), we have 

/ Q\{A, b, c, c)....On(A, b, c, c) e -{SFp{A,b,c,c)+s MFP {U,b,u,^)) 

9=0 J [2?$] e -(SFp(A,b,c,c)+SM F p(U,b,u),u})) 

f VAVbVcVcO x {A,b,c,c)....0 N {A,b,c,c)e- SFP ^ h ~ c ^ 
J VAVbVcVc e -s F p(A,b,c,c) ' 

(22) 

from which one sees that the Faddeev-Popov theory is recovered when 8 = 0. The usual perturbation 
theory is thus re-obtained in the deep ultraviolet region, where the soft parameter 8 becomes negligible. 



{Ox(A,b,c, c)....0 N (A,b,c, c)> 



2.3 Renormalizability and Slavnov- Taylor identities 

2.3.1 Properties of the Gribov propagator and non-renormalization of the soft mass term 

Before showing how the softly broken BRST transformations, eqs. (|18p . can be employed to derive suitable 
Slavnov- Taylor identities, let us give a simple argument in order to convince the reader of the renormal- 
izability of expression (|12p . The reasoning is based on the use of the dimensional regularizatiorjl with 
minimal subtraction and on the following property of the Gribov propagator 

k2 1 294 (23) 



k 4 + 28 4 k 2 & 2 (fc 4 + 20 4 ) ' 

Replacing each Gribov propagator in the one-loop diagrams by expression (|23|) enables us to infer that no 
divergent terms proportional to the soft parameter 8 2 might arise. Let us discuss explicitly the absence 



3 Notice that, in the present case, dimensional regularization in an invariant regularization. 
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of one-loop counterterms of the form 9 Af^A® and 6 2 U^U^. These counterterms would arise, for example, 
from tadpole diagrams in the two point A- A and U-U IP I Green's functions, and would be related to 
the momentum integral 

d«k k 2 (24) 



(2ir) d /c 4 + 20 4 
where d = 4 — e. From identity (|23p . it follows 



r d d k k 2 r d d k 1 r d d k 28* _ 4 r d d k 1 

J (2ir) d k 4 + 28 4 ~ J J2nf k? ~ J {2v) d k 2 (k 4 + 26 4 ) ~ ~ J (27r) d k 2 (k* + 26 4 ) ' ( ' 



where use has been made of dimensional regularization. The right-hand side of equation (J25J) is conver- 
gent in the ultraviolet region by power counting, showing thus that no divergent terms proportional to 8 2 
ariseEl Similarly, one can state that no counterterm of the type # 2 j4^[/^ is required for the one-loop two 
point mixed A-U Green's function, as it follows by noticing that the corresponding Feynman integral is 
ultraviolet convergent by power counting. 

Since no counterterms proportional to 6 2 can arise at one-loop order, it follows that all possible di- 
vergences occurring in the model are those corresponding to 6 2 = 0, in which case the action S reduces 
to that of two completely disjoint expressions, see eq.(|21|). As a consequence, all possible divergences are 
those of the familiar Yang- Mills theory in the Landau gauge, see |X3|, [9] and references therein. Moreover, 
taking into account the mirror symmetry, it follows that the renormalization factors of the mirror fields 
(U^,b a ,H} a ,uj a ) are the same as those of (A,, b a , c a , c a ) which, in the Landau gauge, are given by 19]: 



Zjj = Za 

Zq = Zu = Zc = Z c , (27) 



and 



Z 9 Zfz c = 1 , (28) 

due to the non-renormalization theorem of the ghost-antighost-gluon vertex in the Landau gauge |33j . 
Also, from the absence of one-loop counterterm of the kind 6 2 A'£ l U^, it follows that 

OoA^us, = e 2 Aiu; . (29) 

Therefore, from 

A a — 7 1 ! 2 A a 

US, = Z\ /2 U°, (30) 



we get 



Z 2 = Z~ A X , (31) 



4 It is worth pointing out that the result in eq. (|25p is deeply related to the form of the Gribov propagator. Repeating the 
same argument for a Yukawa type propagator would give rise to a mass dependent divergent term. In fact 

d d k 1 f d d k 1 f d d k A 2 _ 2 f d d k 1 

(2^pfc 2 + A 2 ~ J (2^pF ~ J (27rp k 2 (k 2 + A 2 ) ~" J (27rp k 2 (k 2 + A 2 ) ' ^ > 

which is divergent in the ultraviolet region for d = 4 — e. 
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meaning that the renormalization factoid of the soft parameter 6 2 can be expressed in terms of the 
gluon renormalization factor Zj^. The identity (|23h can be employed to generalize the one-loop argu- 
ment to higher orders, enabling us to argue that the only divergences occurring in our model are those 
corresponding to 6 2 = 0, i.e. those of the Yang-Mills theory in the Landau gauge. 



2.3.2 The Slavnov- Taylor identities 

The existence of a soft breaking of the BRST symmetry does not prevent us to establish a set of Slavnov- 
Taylor identities which are suitable for an all-order algebraic analysis of the renormalizability properties 
of the model [33]. The usual way of proceeding^ is that of introducing an extended action which in- 
corporates all local composite operators entering the soft breaking, by coupling them to a suitable set 
of external sources. The original action is thus recovered when the sources acquire a particular value, 
which we shall refer to as the physical value. The renormalizability of the extended action entails thus 
the renormalizability of the starting action S, eq. (H2l) . 

In order to be able to discuss the non-renormalization property, eq. (|3ip . of the soft parameter 2 , we fol- 
low the procedure outlined in |36] in the case of the study of the composite operator Af^A^ in Yang-Mills 
theory, and introduce a set of sources ( J, rj^, r^, N, M, a®, pV) transforming as 

sJ = , 

st^ = r]^ , 

s-q^ = —d^J , 

spl = o%, 

sal = °> 

sM = N , 

sN = J. (36) 



5 A purely algebraic proof, valid to all orders, of the non-renormalization properties of the soft parameter 2 , eq. (|31[) . is 
under investigation [34] . 

6 We remark here that the non-renormalization of the soft coupling term, eq. (|29[) . seems to be a general feature of theories 
leading to Gribov type propagators, eqs. (|13l) . (|14p . (115p . We observe in fact that a similar result holds in the Gribov-Zwanziger 
theory, for which the corresponding soft coupling term is given by 



7 gf c A; - <p b ; j , (32) 

where 7 is the Gribov parameter and ((f) 1 ^ , 4> b ^) are auxiliary fields, see [9]. The renormalization factors of the Gribov 
parameter and of the auxiliary fields can be found in 35 , and read 

2 r? 2 

7o = Z 7 27 , 

1 be ryl/ 2 1 be 

<^0m = z l /2 <j> b v , (33) 

with 

Z,, = Z-^Z- A 1 ' 4 , < 2 = Z-^zf . (34) 

As a consequence 

2 rabc *a ( ±bc Tbc \ 2 rabc \a ( ±bc ~±bc\ /nr\ 

7offo/ A 0ll _ [ Om - <po^j =1 91 (0 M - M I , (35) 

which expresses the non-renormalization of the soft coupling term of the Gribov-Zwanziger action. 

7 See, for example, the case of the Slavnov- Taylor identities derived in the Gribov-Zwanziger theory |13l |9]. 
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The extended action S turns out to be 



where S 



gext = g 

is defined in eg. (|21|) and Sj is given by 



6»=0 



+Sj 



(37) 



Sj = J (fx (jA«UZ + | J 2 + V ,(c a U« + + r, {s(c a U« + u a A$)\ 



+ 



(fx[N[ 9 r b A c u U a J + gf aco U c u A a y )+M(s( gf aco A c u Uy + gf ac0 U c u A a u u; 



eacbjjc a a , b 



M M 



pacb yt c jjab 



eacbjTC \a , b 



(38) 



where £ is a dimensionless parameter needed to account for the ultraviolet divergences present in the cor- 
relation function ((A°(x)U%(x)) (A h v {y)Ul(y))) . From expression (|38p one easily checks that the starting 
action 5, eq. (fT2|) . is recovered from the extended action S ext when the sources (J, ry^, r M , N, M, a", p°) 
acquire the physical values 



J 



phyt 

(^,T M ,iV,M,<T»,p£ 



so that 



S 



ext 



phys 



phys 

s . 



(39) 
(40) 

(41) 



In order to write down the Slavnov- Taylor identities, a second set of external sources, £1^, L a , L a ), 
properly coupled to the non-linear BRST transformations of the fields (A^, {/", c a , w a ), eqs.(JI8"]). has to 
be introduced [33], namely 



Sn 



(fx -nlD°f(Ay - Q"D?(U)u)° + ^f aoc L a c°c c + ^f aoc L a u co' 



Therefore, it turns out that the action E 

£ = S ext + S n = S 

obeys the Slavnov- Taylor identities 
where 



+Sj + S Q , 



= , 



(42) 

(43) 
(44) 



5(E) 



d 4 x 



SE SE <5£ 5E 5£ SE SE SE <5£ - a SE <5£ 
+ JlT^^" + S^SL*. + + S^ + SQ*~ c W 



+ 



d A x -Juj a 



SE 
Sb G 



5E 
St.. 



SE 



+ N 



SE 



+ J 



SE 



'■5 P * ' At " ' Sr/^ SM " SN 



(45) 



Moreover, besides the Slavnov- Taylor identities (|44jl . the extended action S fulfills other Ward identities, 
which can be found in Appendix [A] These Ward identities turn out to be very helpful for a purely 
algebraic proof [33] of the renormalizability of the model to all orders, a topic which we shall present in 
detail in a forthcoming work [34J. 
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2.4 A few words on the field : relationship with the i-particles 

In order to provide a better understanding of the introduction of the mirror field let us work out the 
relationship between our model, eq. (|12j) . and the i-particles pQ, i.e. the pair of unphysical modes with 
complex masses ±iy/26 2 associated to a confining Gribov type propagator. To this end, let us consider 
the quadratic part of the action S containing the two gauge fields (A®, ?7°), namely 

S quad = J d A x (lA«(-d 2 )A« + ^u;{-d 2 )U; + iV26 2 A°uA , (46) 

where we have already taken into account the Landau gauge conditions, d^Af^ = and d^U^ = 0. 
Expression (|48|) can be cast in diagonal form by introducing the two field combinations 

< = -4 K - u ») ■ ( 47 ) 



V2 

Therefore 



s quad = J d A x ^(-d' + iV^e^ + ^i-d 2 -^ 2 )^ 

which describes in fact the propagation of two unphysical modes with complex masses ±i\^29 2 : 



(48) 



(X^mK-k)) = -Ul{k) + U*(k))(A b v (-k) + U b u {-k))) = 5 a6 fc2+ [5, u - 

WWrfo-k)) = \((A%k) - U^(k))(Al(-k) - U b v {-k))) = S ab k2 _ 1 (Sp, - ffi) .(49) 

These are precisely the i-particles corresponding to the Gribov propagators in eqs. (fl~3j) . (HH), (fT5j) . We 
see thus that the action (|12p has a direct interpretation in terms of i-particles. As pointed out in [lj, 
the advantage of introducing the fields (A^, 77^) relies on the fact that they turn out to be helpful in 
order to construct local composite operators whose one-loop correlation functions exhibit the Kallen- 
Lehmann spectral representation. This feature stems for the observation that the basic momentum 
integral corresponding to the one-loop Feynman diagram containing one propagator of the A-type and 
one propagator of the 77-type, i.e. 

/d^p 1 
(X) 1 {(k - p) 2 + iV20 2 ) ( P 2 -iV2d 2 ) ' (50) 

exhibits a nice spectral representation, as shown in great detail in [I]: 



_ r°° . , , f 1 1 



l(k') - 1(0) = / drp(r) — - - - , (51) 
J2V28 2 \r + k 2 tJ 

where the spectral density 

1 Vr 2 -80 4 

(tr) = 1^2 — — . (52) 

is positive in the range of integratior0. This relevant property enables us to introduce sensible operators 
with good analyticity properties. As an example, let us quote the operator 

O m = {d^K - d v \$ M ~ > ( 53 ) 



8 The subtraction of the factor X(0) in eq. (|5ip is needed to account for the divergent character of expression (|50[1 in four 
dimensions. 
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also extensively investigated in [I], where its two-point correlation function has been shown to be cast in 
the form of a spectral representation with positive spectral functional: 

(0 Xv (k)0 Xv (-k)) = 
p(r) = 

Looking at the expressions of the Gribov propagators in eqs. (fl"3j) . (fl4"J) . (fT5|). we might argue that, in 
the nonperturbative confining infrared region, gluons get affected by nonperturbative long-range effects 
which alter the form of the propagator in such a way that they cannot correspond to excitations of the 
physical spectrum of the theory. Though, these long-range effects have no consequences on the deep 
ultraviolet perturbative region, where the usual interpretation of gluons in terms of partons is reliable. 
We can argue that, in our model, these long-range effects are precisely encoded in the mirror field U^, 
which is only softly coupled to the gauge field A^, through the mixed mass term O^A^U^. This term 
plays a very relevant role in the low energy region, where it gives rise to a confining Gribov type gluon 
propagator. Moreover, in the deep ultraviolet region such a soft coupling is negligible and the two gauge 
fields (A®,U%) belong to two disjoint and non-interacting theories, according to eq. (12ip . In this sense, 
the introduction of the mirror field can be seen as a simple tool to implement a confining mechanism 
for the gluon field A a ^ . Finally, the introduction of the other mirror fields (b a , u a , ui a ) stems from the need 
of preserving the renormalizability of the model. 



2V20 2 r + k 2 ' 



12(N 2 - 1) 



80 4 , 
32vr 2 r 



? 4 + T 2 ) 



(54) 



3 Example of a local composite BRST invariant operator displaying 
the Kallen-Lehmann spectral representation 

Let us face now the construction of a local composite operator whose two-point correlation function 
exhibits the Kallen-Lehmann spectral representation at one-loop order, while being invariant under the 
BRST transformations. To construct such an operator, we observe that the BRST operator s in eq. (|18p 
has non-vanishing cohomology classes |33j which, in the sector of zero ghost number, can be identified 
with BRST invariant local polynomials built with the field strengths F^ U (A), U®(U) and their covariant 
derivatives. Moreover, noticing that the quantity 

F^(x)U^(x) , (55) 

is not BRST invariant^, it follows that the lowest dimensional BRST invariant operators which might be 
introduced and which have well defined transformation properties under the discrete mirror symmetry, 
eq. (|16p . are given by the following two independent combinations 

+ AU {x) = \ {F\x) + U\x)) = \ {F* u (x)F« u (x) + U; v {x)U; v {x)) , (57) 

o- AU {x) = \{f\x)-u\x)) = \{f; v {x)f^{x)-u;^)u;Ax)) ■ (58) 

9 Also here, a suitable subtraction to get rid of ultraviolet divergences is needed, see Q] for details. 
10 An elementary calculation gives 

sF^(x)U^(x) = -gf acb F^ v U^ v c - gf abc F*„U^uj b . (56) 
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As it is easily checked, these two operators have different transformation properties with respect to the 
mirror symmetry (]16p . i.e. 

+ AU (x) -> + <D+ U (x), 

O au {x) -> -0- AU {x). (59) 

Although the operator O au (x) changes sign with respect to the mirror symmetry, its two-point correla- 
tion function (0 AU (k)0 AU (—k)) is, of course, invariant. 

For reasons which will become clear shortly, we shall disregard the operator O au (x), focusing our atten- 
tion to O au (x). We also point out that the two operators O au (x) and O au (x) cannot mix at quantum 
level, as they have different transformation properties under the mirror symmetry, eqs. (|59p . The operator 
O au (x) can only mix with composite operators which are odd under the mirror symmetry. 

The reason why we keep the operator 0~ AU {x) relies on the important property that, to the lowest 
order in the fields, it reduces precisely to the z-particles operator of eq. ([53|) . In fact 

O au (x) = 1 [F« v + U; v ) [F« v - U« v ) = (d n X a u - d v \l) {d^l - d^) + higher order terms . (60) 

As such, we can immediately state that the two-point correlation function {0 AU (k)0 AU (— k)) exhibits 
the Kallen-Lehmann representation at one-loop order. In fact 



poo a/t 2 — R# 4 ('art 4 _i_ T *\ i 
(0- AU (k)0 AU (-k)) = {0 Xn (k)0 Xn (-k)) = 12(iV 2 - 1) / dr^ + j , 2 • 

one— loop j2\/29 2 oZTT T T + K 

(61) 

Concerning the operator O au (x), it turns out that it does not display the same analyticity properties of 
O au (x). This can be checked by expressing it in terms of the fields (A^,^): 

0\u{x) = \{{dX-d^ a af ! + (W-M0 2 )+ higher order terms. (62) 

From expression (i62|) it follows that, at one-loop order, the correlation function (0 AU (x)0 AU (y)) receives 
contributions from terms of the kind ( (d a Xf, — 9^A^)^ (c^-Ap — c^A*)^ ) which gives rise to a Feynman 
diagram with two internal lines of the A^-type. For the resulting momentum integral one gets 



d 4 p p 2 (p — k) 2 + 2(p 2 — pk)' 



(2vr) 4 ((k - p) 2 + iV26 2 ) [f + iy/20 2 ) ' 



(63) 



which does have unphysical cuts along the imaginary axis [Tj. This is a consequence of the fact that, un- 
like expression I(k 2 ) of eq.(|50p. the sign of the complex masses in the denominators of eq. (|63p is the same. 

In summary, we have been able to introduce a BRST invariant local composite operator O au (x), given by 
expression (|58p . whose two point correlation function exhibits nice spectral properties. Other examples 
of higher dimensional BRST invariant composite operators displaying good analyticity properties can be 
easily constructed by inserting an appropriate number of covariant derivatives into expression (|58p . as: 

°dfdu = \ ((Df(A)F^ u ) (D® m (A)F™) - (Df(U)U%) [D<T{U)U™)) . (64) 
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3.1 A short survey on the renormalizability of the operator O 



It remains the hard task of investigating the correlation function (0 AU (x)0 AU (y)) at higher orders, and 
check out if the analyticity properties found at one-loop still hold. Although this analysis is beyond the 
aim of the present work, let us provide here a few preliminary remarks about the renormalization of the 
composite operator O au . 

The property of being BRST invariant will play an important role at the quantum level. Due to the fact 
that the BRST symmetry is softly broken, one expects that the operator 0~2jj will mix with two classes 
of composite operators, as also observed in [11] in the case of the study of the renormalization of the 
composite operator F 2 (x) within the Gribov-Zwanziger theory. 

The first class is given by operators which have the same dimension of O au , i.e. dimension four. These 
operators are BRST invariant, as they are independent from the soft breaking parameter 9 2 . Taking into 
account that 0~ AlJ is odd under the discrete mirror symmetry, eq. ([59|) . and that 0\ v and O au are the 
only two possible non-trivial elements of the cohomology of the BRST operator in the space of the local 
polynomials with dimension four and zero ghost number, it follows that the remaining invariant operators 
with which O au can mix are BRST exact, i.e. they can be written as the s-variation of suitable local 
field polynomials which are odd under the mirror symmetrjjHl- 

The second class of composite operators is given by operators of lower dimension, i.e. of dimension 
two, which depend explicitly on the soft parameter 9 2 . These operators cannot correspond to non-trivial 
BRST invariant elements, as there are no local non-trivial gauge invariant quantities of dimension two. 
Taking into account that O au is odd under the mirror symmetry and that the ghosts and antighosts 
fields (c a ,oj a , c a ,a) a ) cannot enter due to the ghost and antighost Ward identities, eqs. (|103p . (|107p . it 
follows that the only dimension two operator which can potentially show up is 

2 - dm {x) = \ (A«(x)A«(x) - U«(x)U*(x)) . (65) 
It is interesting to observe that this operator has a nice interpretation in terms of i-particles, namely 

°Z-*Jr) = • ( 66 ) 

As such, the correlation function 

{0 2 ^ im {x)0-_ dm {y)) = <(AX)(x) (\U)(y)) , (67) 

will give rise to a momentum integral of the kind of I(k 2 ) in eq. (|5U|) . It will thus display the Kallen- 
Lehmann spectral representation. The same feature holds for the correlation function describing the 
mixing between the two operators O au (x) and ^2-dim( x )- ^ n ^ ac ^> ^° * ne ^ rs ^ or der, 

(02u(x)0-_ dim (y)) = (({d^\ a u - gU£) (« - d uV «)) x (x b A) y ) + higher order terms . (68) 

The Feynman diagram corresponding to expression (|68p will contain one internal line of the A^-type and 
one internal line of the ^-type, giving rise to a momentum integral with a cut structure similar to that of 
I(k 2 ), eq. (|50l) . These properties can be taken as an encouraging evidence that a good analytic structure 

lx For completeness, we mention that composite operators proportional to the equations of motion might also show up. 
However, they do not contribute to the correlation functions [111 112] . 



14 



might survive at higher orders. 



Finally, let us observe that the operator O au {x) can be easily introduced into the extended action 
S, eq. (j43j) . through a BRST invariant source q(x), i.e. 

£ 9 = £ + J d A x q(x)0 AU (x) , (69) 

where the source q(x) is odd under the mirror symmetry, q(x) —> —q(x). The action is seen to obey 
the Slavnov- Taylor identities (|44p 

S(T, q ) = , (70) 

as well as all other Ward identities given in Appendix [A] allowing thus for an all order algebraic analysis 
of the renormalizability of O au (x) along the lines outlined in [T2] . 



4 A possible way to take into account the Gribov copies. The scaling 
type solution for the gluon and ghost propagators 

The action S of our model, eq. (|12p . contains a mass parameter 2 which, till now, has been considered as 
a free parameter. Willing to give a better physical interpretation of the model, a dynamical framework 
for the parameter 2 should be provided, i.e. we should be able to establish a gap equation enabling us 
to express 9 2 as a function of the coupling constant g, as in the case of the Gribov parameter 7 of the 
Gribov- Zwanziger theory [H El [8] . 

A possible way to achieve a meaningful gap equation for 9 2 relies on the observation that the gauge 
fixed action S is plagued by the existence of the Gribov copies [4], as it is apparent from expression (fl~2]) 
where the Landau gauge fixing conditions, d^A® = 0, d^U® = 0, have been employed. 

In [3], it has been suggested that, in order to get rid of the Gribov copies, the domain of integra- 
tion in the functional integral should be restricted to the so-called Gribov region Qa, which is defined as 
the set of field configurations which fulfill the Landau gauge condition and for which the Faddeev-Popov 
operator, M ab (A) = —d^ (5 ab d^ + gf^Afy, is strictly positive, namely 

n A = [A*; d^A" = , M ab (A) = -fy + gf acb A<^ > } . (71) 

The region Qa has been proven to be convex and bounded in all directions in field space [38]. Moreover, 
every gauge orbit passes through Qa [39jlE Its boundary, 8£Ia, where the first vanishing eigenvalue of 
the Faddeev-Popov operator shows up, is known as the Gribov horizon. 

In our case, due to the mirror symmetry (116p . there is also a mirror Gribov region Qjj: 

= {U%; d^l = , M ab (U) = -dp (>d M + gf acb U^j > } . (72) 

12 Nowadays, it is known that the Gribov region Qa is not completely free from Gribov copies [40]. Additional equivalent 
gauge field configurations exist within Qa- To get rid of these additional copies, a further restriction of the domain of 
integration to a smaller region, known as the fundamental modular region, should be implemented. This region is contained 
within the Gribov region Q,a and is known to be free from Gribov copies. Nevertheless, so far, a way to implement the 
restriction of the domain of integration in the functional integral to this region in a local and renormalizable way has not 
yet been achieved. Here, we limit ourselves to the restriction to the Gribov region Ha- 
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Evidently, Q,jj enjoys the same properties of Qa- 



In his original work [3], Gribov was able to implement the restriction to the region Qa by demand- 
ing that the parameter 7 fulfills a particular gap equation, known as the no-pole condition, which follows 
as a consequence of the proper definition of the region £Ia- Here, we shall follow the same route, i.e. we 
shall employ Gribov's no-pole condition in order to establish the gap equation for the parameter 2 . 



Let us proceed by giving a short account on Gribov's no-pole condition, see also [4T] for a pedagogical in- 
troduction. It relies on the observation that, within the regions £7^4, flu, the Faddeev-Popov operators are 
invertible and their inverse, A4~ 1 (A) and .M _1 (C/), are positive definite, as follows from the definition of 
the Gribov regions, eqs. (|71]) . ([72j) . Keeping in mind that M' 1 ^) and M^iU) are nothing but the two- 
point ghost functions, i.e. M.A~ X {x,y) = Gsc(x,y) = (c(x)c(y)), Mu^ix.y) = Gu>w{x,y) = (w(x)w(y)), 
we infer that, within the regions VLa and £lu, Gcc(x,y) and Goui(x,y) stay always positive. 

According to 0], we now parametrize the two point ghost functions Gcc(k), Guu{k) in momentum space 
as 

sab 1 

*=*»=FMV)' (73) 

where the equality Gg(k) = G£t{k) is a consequence of the mirror symmetry f 1 1 6 [) . Gribov's no-pole 
condition amounts to require that the form factor a(k 2 ,0 2 ) is bounded by one, a(k 2 ,8 2 ) < 1, so that 
expression ([73]) cannot have a pole for a non- vanishing value of the momentum k [3]. Expression ([73]) 
stays thus always positive, namely the Gribov horizon 8Qa is never crossed. The only allowed singularity 
is at k = 0, whose meaning is that of approaching the horizon 8Qa, where Gsc(k) is singular due to the 
appearance of zero modes of the Faddeev-Popov operator Ai(A). Of course, the same features hold for 
the mirror two point ghost function G^Lik). 

Following [4], the no-pole condition is implemented by stating that 

a(0,6 2 ) = l, (74) 

which yields the gap equation determining the parameter 6 2 . To the first order, the form factor cr(k 2 , 2 ) 
is easily evaluated and found to be [4j [13] 

a{k ^ } - N 2 -lk 2 J (2<kY (k-q) 2 ^ l9Rl q)) 

2 k fl k u f d d q q 2 ( q^q u 



~ N ° k 2 J (2iry(k-q) 2 (q4 + 20±){ S,MV q 2 J ' (?5) 

where dimensional regularization, d = 4 — e, has been employed. Therefore, at one-loop order, the gap 
equation (i74l) yields 

2 d-1 f d d q 1 



1 = Ng — J WflfTW" (76) 

which enables us to express the parameter 9 2 as a function of the coupling constant g. 

Equation (|76|) also provides us a better understanding of the meaning of the mass paramater 6 2 . We notice 
that, in practice, Gribov's no-pole condition amounts to impose a boundary condition on the two-point 
ghost function G^c(k), i.e. one requires that G^{k) stays always positive, so that the Gribov horizon is 
never crossed. In its original work [4], this boundary condition was implemented by introducing a certain 
massive parameter 7 2 , known as the Gribov parameter, whose value was fixed by the requirement of 
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positivity of the ghost propagator, resulting in a gap equation which has precisely the same form of the 
gap equation ([76]) obeyed by our massive parameter 6 2 , see, for instance, eq.(43) of [I]. This is a relevant 
feature of our model, meaning that no new mass parameters are needed to ensure positivity of the ghost 
propagator, and thus the restriction to the Gribov region Q^. This is achieved by requiring that the 
parameter 9 2 acquires a particular value, determined precisely by the gap equation (j76j) . In other words, 
a unique mass parameter is sufficient in order to ensure positivity of the ghost propagator. As such, our 
massive parameter 9 2 plays the same role of the parameter 7 2 of Gribov's original work, a fact which is 
also apparent by remarking that both parameters give rise to the same gluon propagator, namely 

(A a Jk)A b u (-k)) Gnbov = S^-J^-rfs, k " k " 



(A«(k)A b u (-k)) repllca = ^^^i (<W - ^r) • ( 77 ) 

From the gap equation (|76|) it follows that, for small values of the momentum k, k ~ 0, 

1 - Ck 2 , (78) 



a(k 2 ,0 2 ) 



fc 2 fs0 



for some constant factor d I Thus, for the behavior of the ghost propagators {^(fc), G^ik) near the 
origin in momentum space, we get 



1 



*(*) fe0 = O) fe0 - F . (81) 

showing that imposing the no-pole condition (|76j) yields ghost propagators which are enhanced in the 
infrared. We have thus recovered the so-called scaling solution, i.e. a suppressed gluon propagator which 
vanishes at the origin, eq. (fl3j) . and enhanced ghosts, eq . (|8T]) . This is precisely the type of solution which 
emerges from the Gribov- Zwanziger theory [H El El [8j l4~2l 143] . A scaling type solution is also found in 
the studies of the Schwinger-Dyson equations [21 HSJ 06] , provided the ghost enhancement is required 
by imposing a suitable boundary condition, see [Ml 07J for a recent discussion of this topic. 



5 Clarifying the difference between the Gribov-Zwanziger action and 
the replica model 

5.1 The Gribov-Zwanziger action 

It is worth adding here a few additional remarks in order to clarify in a better way the difference between 
the Gribov-Zwanziger construction and the present model. 



'Equation can be derived by making use of eq. (|76[) , written in the form 



d q 1 fx q^q. 
Therefore, for the form factor a(k 2 , 9 2 ) in eq. (|75[l . we get 



a{k ,6) = 1 + Ng — - / — — - [S^-—\[ - 1 



2 fcfj,fcv 


f d d q 


- 


qp.q v 


k 2 J 




q 2 




f d d q 




q^qv 


k 2 J 


(2n) d 


q 2 



{k - q) 2 J q 4 + 26 4 
1 - N?=& I ^ ( V - ^*-J$ , (80) 



from which eq. (|78|l follows. 
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As it has been already underlined, within the Gribov-Zwanziger approach, a nonperturbative formu- 
lation of Yang-Mills theory is obtained by restricting the domain of integration in the path integral to 
the Gribov region fi^, as expressed by the partition function 

Z GZ = f [DA] 5{dA a ) detid^Dfd^x - y)) €~\ F ^ . (82) 
Jn A 

As shown in [U [6] , the restriction to the region is equivalent to the introduction of the horizon 
function Sh, namely 

Z GZ = J [DA] 5{8A a ) det{d^5\x - y)) e'^^^ F ^ F ^+ S ^) . (83) 

where the horizon ternj 1 "^! Sh is given by 

Sh = 7 4 J d*x h(x) = 7 V J d*xf amc A™ ((— 9 /1 D At )~ 1 ) ad f^A^ (84) 

and 7 2 is the Gribov parameter, which is determined by the so-called no-pole condition [3]. This condition, 
when expressed in terms of the horizon function Sh, takes the form (5j [6] 

(h(x)) = 4(N 2 - 1) . (85) 

At one- loop order, the gap equation (f85|) reads 

2 d-lfd d k 1 
1 = 9 N — J (2nY k* + 2g^ ' (86) 

Taken together, equations (|83p . (|84p . (|85p constitute the so called Gribov-Zwanziger framework for the 
nonperturbative formulation of Yang- Mills theory which takes into account the Gribov issue. 

It is important to emphasize here that the gap equation (|85p is part of the nonperturbative definition 
of the Faddeev-Popov quantization formula within the Gribov-Zwanziger approach. In other words, the 
parameter 7 2 is not a free parameter of the theory. It is generated in a dynamical way, and its expression 
in terms of the gauge coupling constant g is uniquely determined by (|85|) . In particular, one has no more 
the freedom of setting 7 = 0, as one already sees, for example, from the one-loop approximation (|86l) . 
Equations (|83p , (|84|) provide a nonperturbative definition of the quantized Yang-Mills theory only when 
the parameter 7 2 is constrained by the gap equation (J85]). This has a very transparent physical meaning. 
In fact, when 7 2 obeys the gap equation (|85|) it cannot be treated as a free parameter. Instead, it becomes 
a function of the gauge coupling constant g and of Aqcd [4, 5, 6j. As such, expression (1831) has exactly 
the same number of parameters of Yang-Mills theory. This ensures that the restriction of the domain of 
integration in the path integral to the Gribov region £Ia does not bring us outside of Yang- Mills theory. 

5.2 The nonperturbative definition of the replica model 

At this stage, the difference between the Gribov-Zwanziger construction and our approach should be 
manifest. We point out that the nonperturbative definition of the replica model is encoded in the 

14 Although the expression of Sh, eq. (|84[) . is non-local, we remind that it can be cast in local form by means of the 
introduction of a suitable set of localizing fields. The resulting local action can be proven to be renormalizable to all orders 
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following three equations, (jHZD, §B, §B : 



Zrepiica = J [DADbDcDc] [DUDbDuDO] e~ s (87) 
with S being the action of the replica model 

S = Jd A x (\f^(A)F^(A) + \u; u {U)U; u {U) + iV^U^) 

+ j d 4 x (b a d^Al + (fd^{A)c b + b a d^ + uj a dpD ab {U)uj b ^ . (88) 

and 9 2 obeying the gap equation (|74p. namely 

a(0,6< 2 ) = l. (89) 

Taken together, expressions (|87|) . (I§8|) . (|89|) constitute our definition of the nonperturbative version of 
the replica model. In particular, as it happens in the case of the Gribov parameter 7 2 , also the parameter 
9 2 is now completely fixed in terms of the gauge coupling constant g through the gap equation (|89|) . In 
particular, after imposing the gap equation (|74[) . one has no more the freedom to set 2 = 0, as it is 
already apparent from the one- loop approximation (I76p . 

Of course, we are not claiming here that the two partition functions Zqz an d Zrepiica wm gi ye rise 
to the same physics. We remind here that the main motivation to study the replica model is that of 
investigating the consequences of the soft breaking of the BRST symmetry in a confining theory exhibit- 
ing a gluon propagator of the Gribov type. The replica model has to be seen as a useful simpler model 
allowing us to investigate in an analytic way aspects of the gluon confinement. 

We do argue, however, that, taken together, the three expressions (fH7j). ([88]) . ([89]) provide a nonper- 
turbative definition of the replica model which takes into proper account the issue of the Gribov copies, 
as implied by the positivity of the corresponding ghost propagators. As such, the partition function 
Zrepiica does not need to be further restricted to the Gribov region Qa, there is no need for the 
Gribov- Zwanziger construction to be implemented again. This would lead to a complete different model 
which turns out to be plagued by various difficulties, as it will be illustrated in the next section. 



5.3 Combining the replica model and the Gribov-Zwanziger action 

One might argue that, somehow, the replica model and the Gribov-Zwanziger action could be combined 
into a single theory. This would amount to start with two decoupled Yang-Mills theories, i.e. with 6 2 = 0, 
and implement the restriction to the Gribov region fi^ by means of the Gribov-Zwanziger construction. 
At the end, the parameter 6 2 is re-introduced in order to couple in a soft way the two gauge fields A't 
and Up. The final action compatible with the mirror symmetry contains two massive parameters j 2 ,9 2 , 
and reads 

Srepiica-GZ = J #x (^(A)F^{A) + \u^{V)U^{V) + S h (A) + S h (U) + iV2e 2 A«U°\ 

+ [d 4 x(b a dpA^ + c a dpDf{A)c b + b a dpU* + O a dpDfW^ (90) 
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with Sh(A) given in eq.(|84[). Although being a confining model, to our understanding, the action (|90|) 
would represent a completely different model, leading to several difficulties of not easy solution. The first 
one is encoded in the form of the gluon propagator, which takes now a rather complicated form 

<wtc-*» - «- (V - k -t) ■ < 9i > 

This propagator is not of the Gribov type and does not exhibit a simple i-particles decomposition. As a 
consequence, the explicit analytic computation of the spectral densities of the BRST invariant operators 
discussed in the previous sections looks prohibitive. A second difficulty arises at the level of the gap 
equations which should be imposed in order to obtain the values of the two parameters j 2 ,9 2 . Even if 
the no-pole condition 

a(0,9 2 , 7 2 ) = l , (92) 

is required in order to ensure positivity of the ghost propagators, it should be noticed that it will be not 
sufficient to determine both parameters / y 2 ,6 2 - A second gap equation would be required. This is not an 
easy matter. Imposing a gap equation is a highly nontrivial issue, which has to be properly motivated 
by physical and/or geometrical considerations, as in the case of the no-pole condition (|92p which follows 
from the necessity of handling the Gribov copies. At present, it does not look easy to propose a second 
gap equation in order to determine both parameters uniquely. 



6 A refined version of the model. The decoupling type solution 

As already mentioned, the most recent lattice numerical simulations [261 123 US [291 SOI EH [32] point 
towards a gluon propagator which is suppressed in the infrared and which attains a finite non-vanishing 
value at zero momentum, while the ghost propagator turns out to be not enhanced, keeping an almost 

free behavior, G^r(k) ~ -h. 

This behavior, known as the decoupling solution, has also been reported from the analysis of the 
Schwinger-Dyson equations |21 [ 1221 [23] . It can be seen as arising from a different choice of the boundary 
condition for the ghost propagator [16j . 

Within the context of the Gribov- Zwanziger theory, the decoupling solution has been proven to show 
up [201 El EE] when the non-perturbative effects associated to the dynamics of the additional fields 
(<^f, <^ 6 , uj^ , needed to implement the restriction to the Gribov region in a local and renormaliz- 
able way are taken into account. The resulting action is known as the refined Gribov- Zwanziger action 

gang Eg. 

In this section, we intend to show that a refined version of our model can also be introduced. It gives 
rise to a gluon propagator which does not vanish at zero momentum, while preserving locality and renor- 
malizability. This feature relies on the possibility of modifying the starting action 5, eq. (|12|) . by adding 
a kind of Curci- Ferrari mass term |49j . i.e. a mass term whose BRST variation vanishes on-shell, see 
|50 [ I51 [ [52]. In our case, taking into account the mirror symmetry f 1 1 6 1) . this mass term reads 

S m = ^j d A x (A^Al + U;u;) . (93) 
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Accordingly, the action of the refined version of the model is 

Sref = SfP + SmFP + Sq + S m = S + S m 

d'x (±F^(A)F* U (A) + \U« V {U)U; V {U) + iV2e 2 A»U£ + ^ [A-A- + U«U«) 
+ jd*x (b a d^Al + ^d^(A)c b + b a d^U; + Cj a d^{U)j) . (94) 

As it is easily checked, the BRST variation of the mass term S m , eq. (|U3"j) . turns out to be proportional 
to the equations of motion. In fact 

sS m = m 2 f d'x (-A*d fl c a - U«d^) = m 2 j d'x (c J -^f- + , (95) 

so that it vanishes on-shell. This property gives rise to a modified version of the Slavnov- Taylor identities 
which account for the renormalizability of the refined action S re f 

For the gluon propagator, we have now 



,/, , , >v , .,,/, k 2 + m 2 ( k), k„ 



which is of the type already found in the refined Gribov-Zwanziger action [20, 13, 48|^£|. Clearly, expres- 
sion (|96p attains a finite nonvanishing value at zero momentum. 

Even if the issues of a possible dynamical origin of the second mass parameter m 2 in eq. (|93|) and of 
its consequences on the infrared behavior of the ghost propagators are not addressed in the present work, 
the possibility that the decoupling type behavior for the gluon propagator, eq. ([96l) . can be accommodated 
for represents an interesting feature of our model. It is worth observing that, as in the case of the Gribov 
propagator ([1]), the refined propagator (I96D also allows for the introduction of z-particles, namely 

1 (99) 



(k 2 + m 2 ) 2 + 2fl 4 (k 2 + m 2 ) + iy/26 2 (k 2 + m 2 ) - iy/29 2 ' 

In particular, the analysis of the analyticity properties of the BRST invariant operator 0^ v (x), eq. ([58 
generalizes to the case of the refined action (|94p . 



7 Conclusion 

In this work a few aspects related to gluon confinement have been addressed within the Euclidean quan- 
tum field theory. 

15 We mention here that a slightly more general expression for the gluon propagator can be obtained from the generalized 
mass term 

S m , = J d 4 x {^-AlAl + fjU^ , (97) 
which also enjoys the property of being BRST invariant on-shell. For the gluon propagator, we obtain 

= ^ fc4 + (T n 2+M ") fc t: 2 ( 2 , 4 + mV ) (V - k ~f) ■ 08) 
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We have presented a gauge model, eq,(|12|). whose main feature is that of exhibiting a soft breaking 
of the BRST symmetry, encoded in the mass parameter 9 2 , eq. (|lip . This soft breaking is responsible for 
the confining behavior of the propagators of the elementary fields, eqs. (|13j) . (|14p. (|15p . Furthermore, we 
have been able to show that examples of BRST invariant composite operators whose one-loop correlation 
functions exhibit the Kallen-Lehamnn spectral representation can be introduced. 

It is worth to spend a few words on the strategy which has been used in the construction of our model. 
We point out that the mechanism of the soft breaking of the BRST symmetry has played a double role. 
On one side, it is responsible for the confinement of the elementary fields, as expressed by the Gribov 
propagators, eqs. (fT3j) . (fl4j) . (fPoj) . On the other side, it has enabled us to introduce examples of composite 
operators whose one- loop correlation functions turn out to be compatible with the unitarity requirement. 
After that, the existence of the Gribov phenomenon, which is an intrinsic feature of the gauge fixing pro- 
cedure, has been employed to derive the gap equation for the soft parameter associated to the breaking 
of the BRST symmetry. 

To some extent, this framework provides a different way of looking at the Gribov issue. In other words, 
one first pay attention to the construction of a confining model compatible with the requirements of 
renormalizability and of unitarity, encoded in the existence of a set of renormalizable composite opera- 
tors with good analyticity properties. As we have seen, the resulting model turns out to depend on the 
soft parameter responsible for the breaking of the BRST symmetry. As such, it is spoiled of any physical 
interpretation, unless one is able to derive a suitable gap equation allowing for a dynamical determination 
of the aforementioned soft parameter. This is precisely done by invoking the Gribov phenomenon, i.e. 
the existence of the Gribov copies is used as a powerful tool in order to achieve the gap equation for 
the soft parameter. Said otherwise, after one is able to construct a renormalizable and unitary confining 
gauge theory, the Gribov phenomenon is taken as a physical input in order to provide a dynamical frame- 
work for the soft parameter responsible for the breaking of the BRST symmetry, thus making possible 
the contact with QCD. Within this context, it would be interesting to pursue the investigation on the 
analyticity properties displayed by our model, by checking out if these properties can be extended to 
higher orders. 

The main conclusion which emerges from the present work is that the mechanism of the soft break- 
ing of the BRST symmetry might be helpful in order to account for several aspects of a confining theory. 
Certainly, much work is needed in order to achieve a satisfactory understanding of the issue of the BRST 
symmetry versus the Gribov horizon. Nevertheless, the example of the gauge model which we have 
presented here suggests that a possible reconciliation between the BRST symmetry and the presence of 
the Gribov horizon through the mechanism of the soft breaking might provide a useful and interesting 
pathway to gluon confinement. 
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A Ward Identities 



As we have already mentioned in Subsection 12.31 the action E, eq. (jl3j) . fulfills a rather large set of Ward 
identities, which we shall enlist below. Let us first give the full expression of E, namely 

j4„ / ^ jpa I A \ T?<x ( A \ i 1 rm ni\Tja (tt\ , loq <o , ^aa r^ab/ a\J> , Tan jja , -aa r\abi 



d*x l-F^(A)F^(A) + -U* V {U)U* V {U) + b^A* + c^Z^AK + b a d,U* + Q^D^UW 
' d A x (jA%Ul + | J 2 + ^ (c a C7« + u a A«) + r M (J/- 6 ([/»c c c 6 + A^ 6 ) + <?D*(U)u/> + u^^c 6 ) ) 
+ | d 4 x ( 5 / ac6 iV (a^c 6 + ^a; 6 ) - gf acb M ((d,c c )U*c b + A^V + {d^A^J + £/°(9 M c> 6 )) 
+ J d'x (gf abc al (c b U; + JA^ + f (c m c n ^ + W VA£) + <?D^(U) U m + ui b D™(A)c m ^ \ 

+ J (fx (-n*D?(A)<* - n a ^Df(U)u b + 9 -f hc L a c b c c + |/ a6c Z*cA; c ) . (100) 
It turns out that E obeys the following Ward identities: 



the Slavnov- Taylor identities, eq. 

5(E) = 0, (101) 

/" l4 / 5E (5E <5E (5S i5S £E <5E ,„JE 7o <5E T „<5s\ 

C/V"! _ / A^ x I I _i I I ua i_ J „a I 

v ' J [SA^SQ^ 6U£6n° 5c a 5L* 5uj a 5L a <5c a 5u a bb a J 

+ y (-^V + «fe + - + N ui + J m)' (102 > 



the ghost Ward identity [37J E3] 

g m E = A^ , (103) 



where 



d 4 X ( -4- + -r— — + gf aCm C a ^- + 5 / aCm W a X + Ta ( -J— + 



+ / ^ I »/""^ (^ + ^)) • ("*> 



and 



A^ = J d'x {gf acm (^A£ + - LV - + cx^ + " ^ + ^)) • ( 105 ) 

Notice that expression (|105p is linear in the quantum fields. As such, A™ is a classical breaking, 
not affected by quantum corrections [33]. 



• the gauge conditions [33] 

<5S „ .„ aE 
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the antighost equations [33] 



SE 



5H 

5<^ + ^sn?, 



(107) 



the discrete mirror symmetry 



A a - 
A 1 

c a - 
L a - 



J 



of. 



jja 
> b a 

-oo a 

L a 

-»■ J 



M -> M 



w a - 



iV - 

Pi 



A a 
f 1 



L a -> L a 



(108) 



These Ward identities constitute a powerful set up for a purely algebraic investigation [33] of the renor- 
malizability properties of the model to all orders [34] . 
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